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Abstract 

Biochar production and subsequent soil incorporation could provide carbon farming solutions to global climate change and 
escalating food demand. There is evidence that biochar amendment causes fundamental changes in soil nutrient cycles, 
often resulting in marl<ed increases in crop production, particularly in acidic and in infertile soils with low soil organic matter 
contents, although comparable outcomes in temperate soils are variable. We offer insight into the mechanisms underlying 
these findings by focusing attention on the soil nitrogen (N) cycle, specifically on hitherto unmeasured processes of organic 
N cycling in arable soils. We here investigated the impacts of biochar addition on soil organic and inorganic N pools and on 
gross transformation rates of both pools in a biochar field trial on arable land (Chernozem) in Traismauer, Lower Austria. We 
found that biochar increased total soil organic carbon but decreased the extractable organic C pool and soil nitrate. While 
gross rates of organic N transformation processes were reduced by 50-80%, gross N mineralization of organic N was not 
affected. In contrast, biochar promoted soil ammonia-oxidizer populations (bacterial and archaeal nitrifiers) and accelerated 
gross nitrification rates more than two-fold. Our findings indicate a de-coupling of the soil organic and inorganic N cycles, 
with a build-up of organic N, and deceleration of inorganic N release from this pool. The results therefore suggest that 
addition of inorganic fertilizer-N in combination with biochar could compensate for the reduction in organic N 
mineralization, with plants and microbes drawing on fertilizer-N for growth, in turn fuelling the belowground build-up of 
organic N. We conclude that combined addition of biochar with fertilizer-N may increase soil organic N in turn enhancing 
soil carbon sequestration and thereby could play a fundamental role in future soil management strategies. 
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Introduction 

Carbon farming opportunities in temperate agricultural sy.stems 
are limited to no-till farming and conservation agriculture 
strategies. In these systems careful carbon management is required 
to avoid problems of nutrient, particularly N, immobilization by 
the soil microbial community, which can lead to substantial crop 
yield depressions. Biochar application to temperate agricultural 
soils offers the prospect for substantial carbon sequestration 
without the concomitant yield losses or greenhouse gas emissions 
associated with other strategies [1-4]. The pyrolysis derived 
aromatic-macromolecular structure of biochar [5] makes it more 
recalcitrant to microbial decomposition, which means it does not 
immobilize significant quantities of valuable nutrients when added 
to soil, rendering biochar a viable long-term carbon sink option for 
temperate agriculture [6-7]. Promoting novel carbon farming 
options in a wider agronomic context requires facing risk-benefit 
assessments by the soil custodians, namely farmers and other 
environmentally concerned stakeholders. As a result, an in-depth 



understanding of the impact of biochar on the soil N cycle, the 
most economically relevant biogeochemical cycle to farmers in 
mid to high latitudes, will be required to facilitate decision-making 
and allow for improved biochar management and subsequent 
adoption as a carbon sequestering strategy. 

There is increasing evidence suggesting that biochar addition to 
soil impacts a number of processes of the soil N cycle [8] and that 
the impacts of biochar addition are strongly dependent on the 
biochar feedstock and the pyrolysis conditions, in addition to the 
respective soil properties and local environmental and climatic 
conditions [9] . For instance, biochar was shown to decrease nitrate 
leaching [10-11], to increase soil N immobilization [12-13], to 
increase ammonification, the conversion of organic N into 
ammonium [14—15], to increase nitrification, the oxidation of 
ammonium to nitrate [16-17], to promote volatilization, the loss 
of ammonia to the atmosphere [18], and to increase plant N 
uptake [19-21]. The mechanisms posited to explain the differ- 
ences in N cycling processes when biochar is added to soil are 
mosdy direct ones, resulting from the unique properties of the 
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biochar such as its highly porous structure, the large surface area 
and its ion-exchange capacity [22], with associated beneficial 
effects on nutrient retention and water holding capacity. Some N 
cycle processes like NH3 volatihzation can be stimulated by 
inherent biochar properties such as the potential of alkaline 
biochar to increase soil pH [23]. Moreover biochars can have high 
nutrient contents and increase soil nutrient availability directly or 
tiirough the eflFects of priming which may increase the bio- 
availability of soil nutrients [24-25]. In addition, indirect 
mechanisms associated with changes in the soil microbial 
community can be triggered by the properties specific to biochar 
and thereby have strong implications for soil microbial N 
processing [16,26]. To date the majority of biochar studies were 
carried out as short-term pot experiments, many of these focusing 
on poor, acidic and nutrient depleted soils where the addition of 
biochar has the potential to cause largest effects including 
fundamental changes in soil nutrient cycles, often resulting in 
pronounced increases in crop production [27-28]. Long-term field 
trials examining the effects of biochar on intensively managed and 
fertile agricultural soils in temperate regions are therefore urgentiy 
needed to advance our predictive capabilities of biochar - soil N 
cycle interactions [8,20]. 

The impacts of biochar on soil N transformation processes are 
extremely important, given the dominant role N nutrition plays in 
regulating crop production in temperate regions, as it has to be 
replenished either as inorganic N fertilizer or it is cycled back to 
the soil in a complex organic form such as crop residues. In the 
past two decades there has been a major paradigm shift in our 
conceptualization of the soil N cycle [29]. It places greater 
emphasis on the organic turnover processes in the soil N cycle such 
as the depolymerization of organic matter, which is now regarded 
as the rate limiting step of dissolved organic N production and 
subsequently of soil N mineralization and nitrification. As a result 
of the biochar characteristics described above biochar addition to 
soil could have major impacts on the organic N turnover 
component of the soil N cycle. The multitude of simultaneously 
occurring soil N processes (Figure 1) make measuring and 
understanding soil N transformation rates difficult and thus 
requires isotopic tools [30-31]. Using an array of parallel 
estabhshed and novel '^N isotope pool dilution techniques [32- 
33], it should be possible to study each key transformation process 
independendy. These techniques in consortium allow determining 
the impact of biochar on multiple soil processes in a synchronized 
manner and with a high degree of resolution and reliability. In the 
context of biochar addition this issue of complexity is of particular 
relevance as biochar addition could have conflicting impacts on 
different control points of the soil N cycle. 

Our major objective therefore was to investigate the effects of 
biochar on key inorganic and organic N transformation processes, 
on soil C and N pools and on the soil microbial community. A 
hardwood-based biochar (pyrolysis at 500°C for two hours) was 
added in March 20 1 1 to an arable field on calcareous Chernozem 
in Traismauer, Austria. The field trial consisted of four treatments 
replicated four times in circular plots of 6.5 m, comprising NPK 
fertilized plots as controls combined with biochar addition rates of 
0, 24 and 72 t dry mass ha~ . As agricultural systems rely on 
fertilization to grant crop production, fertilized plots were taken as 
controls (NPK, no biochar) and compared to the two biochar 
addition rates i.e. 24 t BC ha~' with NPK amendment (BCIN) 
and 72 t BC ha~' with NPK amendment (BC3N). A treatment 
with 72 t BC ha~' lacking NPK amendment was also setup (BC3) 
to check for nutrient immobilization effects of BC-alone 
treatments. In this study we focused on the control (NPK) and 
the highest biochar treatment (BC3N). Soil samples were taken 
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Figure 1 . Overview of the agricultural soil N cycle. The organic N 
and inorganic N cycles are highlighted. DNRA...dissimilatory nitrate 
reduction to ammonium, DON. ..dissolved organic N, PON... particulate 
organic N such as detritus and soil organic matter, SMB. ..soil microbial 
biomass. Not shown - heterotrophic nitrification which results from 
microbial uptake of e.g. free amino acids, oxidation of reduced N 
compounds and release of nitrate into the soil. 
doi:10.1371/journal.pone.0086388.g001 

four months after biochar application and sampling was repeated 
several times during a period of fourteen months. We applied a 
novel array of isotope pool dilution assays to measure gross rates of 
protein, amino acid, ammonium and nitrate dynamics (Figure 1) 
and molecular techniques to investigate nitrifier (ammonium 
oxidizer) population sizes. In addition we used scanning electron 
microscopy to examine the physical structure of the biochar in 
order to understand its possible role in shaping the soil microbial 
community. 
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Results 

Biochar Porosity 

The cellular structure of wood was well conserved in the biochar 
as seen in Figure S 1 . The mean pore coverage of xylem cross 
sections was 47.6±5.9% (mean±SD) for the biochar applied in 
this field trial. Overall pore size number distribution showed a 
maximum at approximately 0.7 \Lm diameter (related to simple 
xylem pits), and a second and third peak which are closely related 
to bordered xylem pits and parenchyma cells (around 3 |Xm) and 
xylem vessels (15-25 (J.m; Fig. 2, Fig. S2). Beech has a diffuse- 
porous wood structure [34] and early-wood and late-wood had 
distincdy different mean vessel diameters (24±12 and 15 ±9 [Xm, 
respectively), resulting in a decrease in mean pore coverage from 
61.8±6.1% for early-wood to 40.2±4.1% for late-wood. Second- 
ary cortex cells showed a wide range of diameters, with a peak 
ranging from 3 to 5 |J,m (data not shown). On an ar('a basis the 
contribution of small pores decreased and the major contribution 
to total porosity shifted towards larger pores with diameters 
>10 J4,m i.e. xylem vessels (Fig. 2b). 

Effect of Biocliar on General Soil Properties 

The effect of biochar on soil pH and cation exchange capacity 
was minimal. Soil pH decreased by 0.1 units from 7.5 to 7.4 in 
biochar-treated soils (Table 1). Soil cation-exchange capacity 
decreased slightly from an initial 22.5 to 20.8 cmolc kg~' in 
biochar-treated soils but in both cases was dominated by Ca^""" 
(-90%) and Mg^"^ (~7%), with a trend towards an increase of K"^ 
saturation from 1.6 to 4.8% (measured on one composite soil 
sample per treatment, Table 1). 

Bulk Soil C and N Contents and Soil Water Relations 

Biochar amendment resulted in increased soil organic C and 
there was a trend of enhanced total soil N concentration with 
biochar addition. Soil bulk densities and gravimetric soU water 
contents declined while water-filled pore space remained unaltered 
by biochar addition (two-way mixed ANOVA for effects of 
treatment and time, Table 2, Figure 3). Time affected most of 
these and the following soil parameters (Table 2), but emphasis 
here is given on treatment effects and treatment X time 
interactions; the latter might indicate a decreasing or increasing 
biochar-eflfect on soil properties or fluxes. However, all time x 
treatment interactions effects investigated in this study were either 
non-significant or not related to increasing or decreasing biochar- 
efFects over time. Direct time effects were most likely caused by 
seasonal changes, agricultural management and subsequent plant 
growth effects. Soil organic C was significantly higher in biochar- 
treated soils (47.4 mg C g~') than in non-biochar controls 
(18.7 mg C g"'; p = 0.001) (Table 2). Total soil N responded 
analogously, with a slight but non-significant increase due to 
biochar-treatment (1.47 to 1.64 mg N g~'; p = 0.062). Concurrent 
with increases in soil organic C the soil bulk density decreased 
from 1.17 to 1.04 g cm ^ following biochar amendment 
(p<0.001). Along the same lines biochar amendment caused an 
increase in soil porosity from 56.0 to 60.6% (p<0.001) and in 
gravimetric soil water content from 0.140 to 0.163 g H2O g ' dry 
soil (p<0.001). As a result of increased soil water content and soil 
porosity with biochar addition, water-filled pore space (and air- 
fiUed pore space) did not change between controls (29.8%) and 
biochar-treated soils (28.9%, p = 0.229). Interestingly, by compar- 
ing all four treatments we found the biochar-effect to be amount- 
dependent i.e. the magnitude of the effects increased or decreased 
with the amount of biochar added in the field (Table SI). Soil 
organic C increased with biochar rate, bulk density decreased and 



soil porosity and soil water content increased with biochar rate 
(Table SI). 

Extractable C and N Concentrations 

Dissolved organic C in soil extracts was slighdy but significandy 
lower in the biochar treatment compared to the controls (111 vs. 
123 |Ig C g"' respectively; p = 0.010) (Figure 3, Table 2). 
Concentrations of extractable ammonium were low compared to 

extractable nitrate and not significantly different between the 
control and biochar treatments, with 1.66 and 1.55 |Xg NH4''^-N 
g ' respectively (p = 0.383). In contrast soil nitrate concentrations 
in the biochar (BC3N) treated soils were significantly lower than 
those from the control soils (p = 0.001), with 8.40 and 14.4 |a,g 
NOj^N g~' respectively. Dissolved organic N concentrations in 
the soil extracts decreased, though not significantly, from 9.62 in 
bio[:har-treated soils to 8.35 |J.g N g in txjntrol treatments 
(p = 0.187). The ratio of dissolved organic C: dissolved organic N 
varied around 16:1 and was not affected by biochar addition 
(p = 0.909). Nitrate dominated the extractable N pool, accounting 
for between 67% (BC3N) and 76% (controls). The decrease in soil 
nitrate therefore translated into a significant decline of total 
dissolved N from 25.0 to 17.0 |J,g N g~ in the biochar treatment 
(p = 0.001). Free amino acids were low in concentration (around 
0.7 (Xg N g~') and were not affected by treatment (p = 0.110). 
Measured only once (July 201 1, four months after the start of the 
field experiment) we found no significant differences in the 
concentrations of different size classes of dissolved organic N 
between the two treatments, with extractable high-molecular 
weight organic N concentrations around 3.3 |a,g N g~' and low- 
molecular weight organic N concentrations of around 6 Hg N g ' 
(Table 2). Free amino acids accounted for approximately 12% of 
the low molecular weight organic N pool and comprised 8% of 
dissolved organic N in both treatments. Around 40%) of the 
dissolved organic N and all the high-molecular weight organic N 
fraction (>3 kDa; 105-107'X>) was protein, although this mea- 
surement is based only on the July 20 1 1 sampUng. The negative 
effect of biochar on concentrations of dissolved organic C, total 
dissolved N and nitrate was consistent even when all four 
treatments were analyzed (Table SI), with decreases in dissolved 
organic C, nitrate and free amino acids with increasing biochar 
addition rate (Table SI). 

Soil N Transformation Rates 

The fuU set of gross N transformation processes was measured 
once, and one-way ANOVA results are presented in Table 3 and 
Figure 4; nitrate fluxes were measured on three occasions 
throughout the seasons (Table 3, Figure 4, see circles). Measured 
rates of protein influx into the extractable organic N pool were 
significandy (p = 0.016) higher in the NPK treatment (80.0 [Ig N 
g~' d~') and almost double the rates in the biochar treatment 
(44.2 Hg N g"'d''). Biochar showed a trend (p = 0.097) of 
depressing protein efflux rates out of the extractable organic N 
pool, with rates in the biochar treatment (53 ft^ N ^~ ' d~ ') being 
half that of the control treatment (1 15 |Xg N g d ). The decline 
in amino acid transformation rates was even more pronounced. 
Amino acid influx rates into the organic N pool were significantly 
(p = 0.0004) lower in biochar treated soils (5.4 |J.g N g ' d '), 
more than fivefold lower than those in the control treatment 
(28.5 |i,g N g ' d ') and so too were amino acid consumption 
rates. Amino acid consumption rates declined more than six-fold 
from 50.3 |J,g N g ' d 'in the control treatment to 7.9 |J,g N 
g ' d 'in the biochar treatment (p = 0.039). Gross mineralization 
rates were not significandy different between the two treatments 
(p = 0.616) nor were ammonium consumption rates (p = 0.526; 
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Figure 2. Pore size distribution of beech wood biochar as determined by scanning electron microscopy. In total, mean diameters of 
7715 pores were measured in 6 longitudinal and 7 cross sections of charcoal pieces, each image covering 0.35 mm^. Pore size distributions are 
presented (1, top panel) as counts mm^'^ per 1 |im size class (top) and (2, middle panel) on an area basis, i.e. |a.m^ mm"'^ per 1 )im size class. Typical 
mean diameters (soil bacteria, soil fungi) and diameter ranges (arbuscular mycorrhizal fungi; protozoan grazers such as flagellates, amoebae, ciliates; 
and mesofaunal grazers and predators such as nematodes, mites and collembola) are presented color-coded (bottom panel) and with range lines in 
the middle panel. 

doi:1 0.1 371 /journal.pone.0086388.g002 



Figure 3), and both processes ranged between 1.5 and 2.1 Hg N 
g ' d In contrast, gross nitrification and gross nitrate 
immobilization were strongly enhanced by biochar application. 
Gross nitrification rates in the biochar treatment (10.2 pg N 
g~ d~ ) were fourfold higher than the rates in the control 
treatment (2.4 pg N g d ), and significantly different 
(p<0.001). A similar trend (p = 0.0573) was evident for nitrate 



consumption rates which almost doubled from 6.8 to 11.3 pg N 
g~' d~' in the biochar treatment. 

Rates of nitrification and nitrate consumption were re-measured 
twice. Over the course of time nitrification rates were consistently 
higher (threefold) in the biochar treated compared to the control 
plots (j3<0.001, Table 3). The same pattern was apparent for 
nitrate consumption rates which were consistently negatively 
affected by biochar amendment (p = 0.001). On the sampling date 
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Table 1. Soil properties in control (NPK) and biochar-treated 
(BC3N) soils of the Traismauer biochar field trial. 



Parameter 


Unit 


Control (NPK) 


Biochar (BC3N) 


pH (CaCb) 




7.5 


7.4 


CaCOj 


% 


15.8 


15.2 


Humus 


% 


2.4 


18.1 


Total N 


% 


0.148 


0.203 


P (CAL) 


mg kg"^ 


49 


84 


Ptot {3cid digest) 


g kg-^ 


5.46 


5.54 


Sand 


% 


18.3 


n.d.' 


Silt 


% 


57.2 


n.d. 


Clay 


% 


24.5 


n.d. 


CEC 


cmol kg"^ 


22.5 


20.8 


Ca (CEC) 


cmol kg"^ 


20.7 


18.2 


Mg (CEC) 


cmol kg"^ 


1.46 


1.53 


K (CEC) 


cmol kg"^ 


0.36 


0.99 


Na (CEC) 


cmol kg"^ 


<C.C4 


<0.04 


Al (CEC) 


cmol kg"^ 


<0.06 


<0.06 


Fe (CEC) 


cmol kg"^ 


<C.C1 


<0.01 


Mn (CEC) 


cmol kg"^ 


<0.01 


<0.01 


H (CEC) 


cmol kg"^ 


0.002 


0.002 


Fe (EDTA) 


mg kg"^ 


40 


67 


Mn (EDTA) 


mg kg"^ 


107 


128 


Cu (EDTA) 


mg kg"^ 


7.2 


7.1 


Zn (EDTA) 


mg kg"^ 


2.3 


7.5 



All parameters were measured on composite samples collected in July 2011 

from the two treatments (n = 1). n.d. - not determined. 

^Soil texture is not shown for biochar treated soils. At high contents humus 

biases the particle size distribution and is therefore conventionally decomposed 

by prolonged H2O2 (or similar oxidant) treatment prior to soil texture analysis; 

biochar/charcoal is largely resistant against this oxidant treatment. 

doi:l 0.1 371 /journal.pone.0086388.t001 

where we had gross nitrification rates measured for all four 
treatments nitrification rates were significantly higher in all the 
biochar treatments (p<0.001, Table SI) with rates increasing 
significantly with increasing rates of biochar application. 

Microbial Connmunity Size and Abundance of Archaeal 
and Bacterial Ammonia-oxidizers 

Soil DNA concentration was significantly higher in the biochar 
treatment (18.0 |a,g DNA g~' soil) compared to the controls 
(12.4 Hg DNA g"'; p = 0.042) (Figure 5, Table 2). Ammonia- 
oxidizing archaea (AOA) outnumbered bacterial ammonia- 
oxidizers (AOB) by a factor of 3.2 to 3.6, in both treatments 
(p = 0.402). At the same time archaeal amoA gene copy numbers 
(AOA) and bacterial amoA gene copy numbers (AOB) expressed 
on a dry soil mass basis increased through biochar application 
(AOA: P = 0.045, AOB: p = 0.020). AOA increased 1.5-fold from 
1.68x10' to 2.47x10' amoA copy numbers g ' soil and AOB 
1.7-fold from 5.01 x 10** to 8.45x10** amoA copy numbers g ' soil 
through biochar amendment. As a result of increasing soil DNA 
concentration and amoA copy numbers per gram soil we did not 
find a significant effect of biochar on the abundance of AOA 
(p = 0.637) and AOB (p = 0.436) expressed on a DNA basis, i.e. the 
relative share of ammonia oxidizers in the microbial community 
did not change according to treatment. Across all four treatments 



we did not find an effect of biochar rate on either soil DNA 
concentration or any measure of AOA or AOB abundance (Table 
SI). However, archaeal and bacterial amoA gene copy numbers 
were strongly positively related (R^ = 0.83, p = 0.01 14; data not 
shown). Moreover, we found a trend towards a positive 
relationship between both, AOA or AOB numbers on a dry soU 
mass basis and gross nitrification rates (AOA: = 0.65, p = 0.051, 
AOB: R^ = 0.60, p = 0.072; Figure 6). The "weakness" of this 
relationship derived mainly from the control treatment where 
nitrification rates were relatively constant between the two 
measurement dates while amoA copy numbers varied more 
strongly. Nitrification rates were significantly positively correlated 
with the total number of amoA copy numbers on a soil dry mass 
basis (AOA plus AOB; R^ = 0.66, p = 0.049, data not shown). 

Discussion 

Soil Organic N Cycling 

Biochar addition had a strong negative effect on soil organic N 
cycling rates as it reduced protein production and protein 
consumption and concurrently decreased free amino acid 
production and consumption (Figure 4). This is important given 
that protein depolymerization plays the predominant role in 
organic N recycling with approximately 50% of soil organic matter 
and >90% of plant detritus being composed of protein [35-36], 
and with proteins representing around 40% of dissolved organic N 
in this study. Proteins are therefore considered the most important 
substrate for microbial N mobilization and N use [33], and thus 
any change in their dynamics strongly impacts whole-system N 
cycling. The influx into the KjSO^ extractable protein pool 
(representing a labile part of the total protein pool in soils) is 
mainly fuelled by microbial turnover and protein desorption 
processes while efflux from this pool relates to protein depolymer- 
ization or protein adsorption (Figure l).The concomitant reduc- 
tions in free amino acid turnover were most likely a direct 
consequence of the reduced protein-N cycling in biochar amended 
soils. Amino acid production is driven by depolymerization of 
accessible proteins through extracellular proteases and peptidases 
which release oligopeptides and free amino acids that are then 
taken up by microbes (free amino acid consumption) [33]. Though 
oligopeptides also represent a significant source of N for microbes 
and plants [37] we did not measure their dynamics, which would 
have been expected to change accordingly. We offer two possible 
explanations for the observed negative effects of biochar on 
organic N (protein) dynamics related to the surface properties of 
biochar and its unique pore structure. 

(1) The observed decline in protein production rates might result 
from protein adsorption to biochar surfaces. The reactive 
surface of biochar has been shown to adsorb numerous low 
and high molecular weight organic molecules [38], including 
proteins and enzymes [39]. Hence, the reduced motility and 
accessibility of proteins in the biochar treated soils could 
explain the decrease in protein influx rates observed with 
resultant implications for amino acid dynamics. Very strong 
H-bonding between proteins and polyphenols are well known 
[40], and the abundant phenolic groups on low temperature 
biochars [41] provide ample ligands for strong sorption of 
proteins and peptides [9] . 

(2) Biochar could also provide structural protection, by the 
entrapment of organic matter within its inherent pore 
network. The "mesopore protection hypothesis" by Zimmer- 
man et al. [42] suggests that small organic molecules can be 
occluded within mineral mesopores (2-50 nm diameter). 
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Table 2. Results of ANOVA of soil physicochemical data for the factors treatment and time in control (NPK) and biochar (BC3N) 
treatments. 



Pools 


Unit 




NPK 


BC3N 


Treatment 




Time 




Interaction 










mean 


mean 


F 


P 


F 


P 


F 


P 


^org 


mg C g"' DW 




18.66 


47.43 


(1, 6) 43.7 


0.001 


(1, 6) 0.0 


0.951 


2.5 


0.165 


Nto, 


mg N g"' DW 




1.47 


1.64 


(1, 6) 5.3 


0.062 


(2, 12) 37.5 


<0.001 


0.3 


0.778 


BD 


g cm"^ 




1.17 


1.04 


(1, 6) 48.4 


<0.001 


(1, 6) 55.3 


< 0.001 


0.3 


0.595 


Grav. WC 


g H2O g"' DW 




0.14 


0.17 


(1, 6) 46.7 


<0.001 


(6, 36) 218.7 


<0.001 


6.2 


< 0.001 


Porosity 


% 




56.03 


60.61 


(1, 6) 48.6 


<0.001 


(1, 6) 55.3 


<0.001 


0.3 


0.594 


WFPS 


% 




29.79 


28.92 


(1, 6) 1.8 


0.229 


(6, 36) 215.4 


<0.001 


2.3 


0.059 


APRS 


% 




70.21 


71.08 


(1, 6) 1.8 


0.229 


(6, 36) 215.4 


< 0.001 


2.3 


0.059 


DOC 


Hg C g"' DW 




122.56 


111.21 


(1, 6) 13.6 


0.010 


(5, 30) 205 


<0.001 


2.9 


0.028 


DOC:DON 






16.10 


15.92 


(1, 6) 0.0 


0.909 


(1.6, 9.9) 14.4 


0.002 


0.8 


0.454 


DON 


|ig N g"' DW 




9.62 


8.35 


(1, 6) 2.2 


0.187 


(2.3, 13.9) 28.5 


<0.001 


0.3 


0.782 


NO3" 


Hg N g"' DW 




14.43 


8.40 


(1, 6) 39.6 


0.001 


(6, 36) 95.9 


< 0.001 


8.6 


< 0.001 


NH4+ 


[ig N g"' DW 




1.66 


1.55 


(1, 6) 0.9 


0.383 


(2.6, 15.5) 34.8 


<0.001 


1.0 


0.417 


FAA 


[ig N g"' DW 




0.77 


0.70 


(1, 6) 3.5 


0.110 


(1.6, 9.8) 3.2 


0.091 


0.9 


0.429 


TDN 


Hg N g"' DW 




24.95 


17.03 


(1, 6) 33.3 


0.001 


(5, 30) 97.2 


<0.001 


1.3 


0.285 


HMW-N„,g 


[ig N g"' DW 




3.24 


3.47 


(1, 23) 1.1 


0.306 










LMW-N„,g 


Hg N g"' DW 




6.78 


5.68 


(1, 23) 0.7 


0.401 










Protein 


[ig N g"' DW 




3.46 


3.65 


(1, 23) 1.4 


0.253 










DNA 


Hg g"' DW 




12.40 


18.04 


(1, 6) 6.6 


0.042 


(1.1, 6.9) 0.1 


0.774 


2.6 


0.154 


AOA/DNA 


amoA copy no ^g" 


' DNA 


1.41E6 


1.36E6 


(1, 6) 0.2 


0.637 


(2, 12) 8.2 


0.006 


3.3 


0.070 


AOA/soil 


amoA copy no g"^ 


soil DW 


1 .68E7 


2.47E7 


(1, 6) 6.3 


0.045 


(1.1, 6.6) 3.5 


0.103 


0.9 


0.378 


AOB/DNA 


amoA copy no ^g" 


' DNA 


419075 


461583 


(1, 6) 0.7 


0.436 


(1.2, 7) 7.2 


0.028 


1.6 


0.246 


AOB/soil 


amoA copy no g"^ 


soil DW 


5.01 E6 


8.45 E6 


(1, 6) 9.9 


0.020 


(2, 12) 4.5 


0.035 


1.3 


0.319 


AOA:AOB 






3.64 


3.19 


(1, 6) 0.8 


0.402 


(2, 12) 0.4 


0.681 


0.0 


0.970 



Data were analyzed by two-way mixed ANOVA. Data for protein, HMW-Norg and LMW-Norg were available only for one sampling date and therefore were analyzed by 
one-way ANOVA. Abbreviations: Corg, soil organic C; Ntot. total soil N; BD, soil bulk density; Grav. WC, gravimetric soil water content; WFPS, water filled pore space; AFPS, 
aire filled pore space; DOC, dissolved organic C; DON, dissolved organic N; NO3", nitrate; NH4\ ammonium; FAA, total free amino acids; TDN, total dissolved N; HMW- 
Norg/ high molecular weight organic N; LMW-Norg, low molecular weight organic N; DNA, soil DNA content; AOA/DNA, archaeal amoA copy numbers on DNA basis; AOA/ 
soil, archaeal amoA copy numbers on dry soil basis; AOB/DNA, bacterial amoA copy numbers on DNA basis; AOB/soil, bacterial amoA copy numbers on dry soil basis. 
doi:l 0.1 371 /journal.pone.0086388.t002 



physically protecting organic molecules such as proteins 
(<10nm) and amino acids (<2 nm) from breakdown by 
excluding decomposers and their extracellular enzymes. 
Mesopore surface area and volume are mainly determined 
by gas sorption analysis but are not targeted by scanning 
electron microscopy, and therefore have not been measured 
in this study. Nevertheless, low-temperature wood biochars 
produced at <650°C, just as the Romchar used in this study, 
have been shown to have a low mesopore volume and surface 
area as measured by gas sorption [39] . Physical protection of 
organic matter in mesopores by occlusion from decomposing 
microbes may thus have played a minor role in this study. 

Just like the decrease in protein production, the decline in free 
amino acid production in biochar treated soils could also be 
triggered by protease and peptidase adsorption to or occlusion in 
biochar. Enzyme adsorption to solid support materials (e.g. soil 
matrix or biochar) has been reported to decrease soil enzymatic 
activity [43], and adsorbed substrates such as proteins and 
oligopeptides are also less accessible and bioavaHable for further 
breakdown to free amino acids. Decreased aminopeptidase activity 
has been observed in experiments with fast pyrolysis biochar of 
switchgrass and was ascribed to sorption reactions between 



substrate and biochar [44]. Another study found accelerated N- 
mineralization of the recalcitrant soil organic N pool in arable soils 
treated with maize biochar [14]. The inconsistent effects of 
biochar on soil enzymes may be the result of several factors such as 
biochar feedstock, pyrolysis conditions, soil type as well as the 
respective enzyme and/or substrate [44]. For example, contrasting 
results obtained in our experiment and that of Nelissen et al. 
(2012) [14] could be explained by structural differences between 
the biochars used in these two studies, leading to variable enzyme 
diffusion impedance. Indeed, the Romchar employed here 
maintained its native wood-like structure according to the 
scanning electron microscopy images (Figure SI) while non-woody 
materials such as maize biochar have an amorphous structure with 
litde evidence of retained pores [45] . 

Soil N Mineralization 

Given that biochar application significandy decelerated the 
production and consumption of low molecular weight organic N in 
our field experiment, we would have expected a decrease in soil N 
mineralization rates. Somewhat surprisingly biochar had no 
significant effect on gross N mineralization and gross ammonium 
immobilization rates in our study, in line with a laboratory 
incubation of cultivated black Chernozem soil with wheat biochar 
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Figure 3. Soil C and N pool sizes for control (NPK) and biochar 
(BC3N) treatments. Bars represent means ±1 SE (n = 4) from seven 
measurements between July 201 1 and September 2012 for extractable 
compounds, and from two measurements (September 2011 and 2012) 
for soil organic C and total soil N. Open bars, control treatment (NPK); 
grey bars, biochar treatment (BC3N). Units in mg C g"^ dry soil and |ig 
N g"^ dry soil. P values are from two-way mixed ANOVA (see Table 2). 
doi:10.1371/journal.pone.0086388.g003 

[46], but in contrast to observed accelerated rates of gross N 
mineralization of the recalcitrant soil organic N pool in an arable 
soil treated with maize biochar [14]. Interestingly mineralization 



(respiratory use) of amino acids and peptides, measured as the 
uptake of '''C-labelled organic N and subsequent '*C02 release, 
was not affected by application of a wood-based biochar [15]. 
Divergent, positive, neutral or negative effects of biochar on soil N 
mineralization may be the upshot of different biochar types along 
with impacts on the C and N status of the soil microbial 
community. Carbon-rich biochar, poor in nitrogen, such as beech 
wood biochar (C:N ratio ~200:1) might promote N limitation, 
leading to the retention of catabolically produced ammonium in 
the N-limited microflora, which therefore results in a decrease of 
the amount of ammonium released to the soil [29]. Conversely, 
nitrogen-rich biochars with low C:N ratios such as maize-biochars 
[14] or manure-biochars [47] promote microbial C limitation, 
causing the excess of N to be mineralized and therefore N 
mineralization rates to increase. In our case the ammonium 
released (N mineralization) relative to organic N taken up by 
microbes (free amino acid consumption) was higher in the biochar 
treated soils than in controls (Figure 4), indicating that biochar 
addition did not lead to microbial N limitation as a consequence of 
the increased organic C content in soil. This suggests that the 
majority of organic C in the biochar was not metabolized or 
assimilated by the soil microbial community, and confirms that 
biochar is an efiicient tool for long-term C sequestration in soil, 
with half-lifes >100 years [48]. The absence of microbial N 
limitation in the biochar treatment might be directly related to the 
field applications of fertilizer and might also be related to the 
slightly lower concentrations of dissolved organic C in the biochar 
treated plots. 

Soil Nitrification 

Apart from the negative elfect that biochar had on soil organic 
N cycling, the most striking result was a strong and persistent (1 .5 
years) stimulation of gross nitrification. Several studies have 
examined the impact of biochar on soil net nitrification [8] but 
only few have reported gross nitrification rates, which allow the 
simultaneous analysis of both, productive and consumptive 
processes for nitrate. Biochar addition increased gross nitrification 
rates in two forest soils and one arable soil [14,26,49] but had no 
effect in other experiments with grassland and arable soils [46,49] . 
Several hypotheses have been raised to explain biochar-related 
(positive) effects on nitrification rates. 

One of the dominant theories to explain the increase in 
nitrification rates following biochar amendment is the liming effect 
[14,26]. Biochar commonly has a neutral to alkaline pH [50] and 
therefore may increase the pH of acidic soils with a concurrent 
positive response of soil nitrifiers whose pH optimum is sUghtiy 
acidic to neutral pH [51-52]. However the pH of the soil used in 
this study was slightly alkaline and did not increase but rather 
decreased in response to the addition of biochar. 

Equally we can disregard the tenet that biochar addition leads 
to an increase in nitrification due to adsorption of nitrification 
inhibitors. In contrast to forest and grassland soils [49,53-54], 
intensively managed agricultural soils usually contain low amounts 
of naturally occurring nitrification inhibitors [55] and their effects 
might have been of minor importance in this experiment. 

Nitrification rates have been shown to be strongly controlled by 
the ammonium supplied by soil N mineralization in a wide variety 
of ecosystem types including woodland, grassland and agricultural 
soils [56] . In this study, biochar amendment neither had an effect 
on gross ammonium production rates nor on the pool size of 
extractable ammonium, which obviates a simple substrate-driven 
explanation of increased nitrification rates in this field experiment. 
Whereas nitrification rates in the control soils were of the same 
order of magnitude as the mineralization rates, biochar treated 
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Table 3. Results of ANOVA of soil N transformation rates for the factor treatment in control (NPK) and biochar (BC3N) treatments. 



Processes 


Unit 




NPK 


BC3N 


Treatment 




Time 




interaction 








mean 


mean 


F 




P 


F 


P 


F P 


One date 


Prot_GP 


mg N l<g" 


d"' 


79.99 


44.15 


(1, 


1 1 ) 8.4 


0.016 








Prot_GC 


mg N l<g" 


d"' 


1 1 5.06 


52.37 


(1, 


1 1 ) 3.4 


0.097 








AA_GP 


mg N l<g" 


d"' 


28.54 


5.36 


(1, 


1 1 ) 26.7 


0.000 








AA_GC 


mg N l<g" 


d"' 


50.25 


7.91 


(1, 


11) 5.6 


0.039 








N min 


mg N l<g" 


d"' 


1.88 


2.10 


(1, 


1 1 ) 0.3 


0.616 








NH4+_GC 


mg N l<g" 


d"' 


1.86 


1.52 


(1, 


1 1 ) 0.4 


0.526 








Nitrification 


mg N l<g" 


d"' 


2.42 


10.16 


(1, 


1 1 ) 50.6 


<0.001 








NOj^.GC 


mg N l<g" 


d"' 


6.83 


11.27 


(1, 


1 1 ) 4.6 


0.057 








Three dates 


Nitrification 


mg N l<g 


d-' 


3.38 


9.09 


(1, 


6) 82.2 


<0.001 


(2, 12) 0.4 


0.698 


3.6 0.061 


NOj^.GC 


mg N l<g" 


d-' 


4.92 


11.01 


(1, 


6) 35.2 


0.001 


(2, 12) 1.7 


0.225 


0.6 0.571 



All transformation rates were measured and analyzed by one-way ANOVA for July 201 1 data ("one date"). Nitrification and nitrate consumption was measured three 
times for NPK and BC3N treatments ("three dates") and analyzed by two-way mixed ANOVA. Abbreviations: GP, gross production (Influx); GC, gross consumption 
(efflux); Prot, protein; AA, amino acid; N mIn, N mineralization; NH4^, ammonium; NO3", nitrate. 
doi:1 0.1 371 /journal.pone.0086388.t003 



soils exhibited nitrification rates four to five times higher than 
mineralization rates. This dichotomy between ammonium supply 
and nitrification rates could point to increased importance of 
heterotrophic nitrification fuelled by organic N compounds in 
biochar-treated soils. However heterotrophic nitrification mainly 
occurs in grassland and forest soils and especially in acidic soils 
[57-58], although this process has been recently detected in acidic 
arable and neutral soils [59-60]. This nitrification pathway will 
have to be investigated in follow-up incubation experiments, using 
acetylene as an inhibitor of autotrophic (archaeal and bacterial) 
nitrification [61] which does not affect heterotrophic nitrification 
[62-63]. 

Increased nitrification rates in soils treated with biochar could 
be attributed to the alteration of the soil water status, including its 
distribution and associated changes in soil oxygenation, since 
nitrification is dependent on soil O2 availability. The soil air/water 
balance is essential for nitrification, which has been found to peak 
at 60% water-fdled pore space [64]. Altered soil moisture 
conditions have been reported in various biochar studies [17,65- 
66]. However, in this study air-fdled pore space was highly 
variable and on average greater than 70%, independent of 
treatments; therefore higher nitrification rates (Figure 4) could not 
be explained by improved oxygen availability. 

While none of the proposed drivers of soil nitrification could 
explain the observed pattern of enhanced nitrification in biochar 
treated soils, the size of the ammonia-oxidizing community 
increased with biochar amendment (Fig. 5). Recently, the 
oxidation of ammonia to nitrite, the first and rate limiting step 
of nitrification, was shown to be performed by prokaryotes in both, 
the bacterial and archaeal domains [67]. Both ammonia oxidizing 
archaea (AOA) and ammonia oxidizing bacteria (AOB) abun- 
dances were significantly higher in the biochar treated soils, AOA 
outnumbering AOB in both treatments. There are numerous 
studies showing that AOA are often significantly more abundant 
than AOB in soils of different origins [68-69]. However the 
relative importance and contribution of these different groups to 
the nitrification process in natural environments is still unclear and 
to our knowledge there are no published studies available on the 



interaction of biochar and AOA, at least in temperate ecosystems. 
AOB abundances have been found to increase after natural forest 
fire events most likely as a consequence of the liming effect of 
charcoal and its sorption of nitrification inhibitors [26] . In a pot 
experiment the AOB community structure was unaffected by the 
sole addition of Eucalyptus wood biochar but changed when 
biochar was added together with an inorganic or organic N source 
[70], a possible result of substrate-limited nitrification as the 
ammonium concentrations declined in sole biochar treatments. 
The stimulation of both, AOA and AOB abundances through 
biochar addition in this field experiment was only significant when 
amoA gene copy numbers were expressed on a soil mass basis but 
not when expressed on a soil DNA basis. The microbial 
community size, indicated by total soil DNA content, increased 
with biochar addition revealing a general positive effect of biochar 
on the whole microbial community rather than a specific 
stimulation of the autotrophic nitrifiers. This unspecific stimulation 
of microbes and concomitant increase of nitrifiers is corroborated 
by the lack of biochar effects on more specific soil nitrification 
controls such as soil pH, ammonia and oxygen availability, and 
presence of nitrification inhibitors. The biochar effect was also 
unspecific in terms of stimulation of AOA and AOB; both AOA 
and AOB abundances increased with biochar and were strongly 
positively related across all soil samples independent of treatment. 
In addition ratios of AOA: AOB were unaffected by the biochar 
amendment. 

AOA and AOB abundances were broadly (P<0.1) correlated 
with gross nitrification rates, suggesting that bacterial and archaeal 
abundance increased with increasing biochar addition and 
providing a functional link to the process data (Fig. 6). This and 
the total DNA data add weight to an underlying hypothesis, i.e. 
that the structural nature of the biochar provides a protective 
habitat for bacterial and archaeal growth. It is conceivable, that 
biochar promotes bacterial and archaeal community size by 
enhancing the colonizable surfaces in the soil and providing 
habitable protected pores. The unique pore-size distribution 
typically found in biochar with pores ranging from 1 nm to 
100 |a,m [71] is believed to provide various beneficial functions to 
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Figure 4. Soil N transformation rates for control (NPK) and 
biochar (BC3N) treatments. Bars represent means ±1 SE (n = 4) of 

measurements from July 201 1 . Open bars, control treatment (NPK); grey 
bars, biochar treatment (BC3N). Open circles represent means ±1 SE 
(n=4) of nitrification rates from three measurement dates. Left panels 
represent gross influx rates into the target pools (production rates), 
right panels gross efflux rates (consumption rates). Units in [ig N g"^ 
dry soil d"\ equivalent to mg N kg"' d"\ P values are from one-way 
ANOVA (all rates) or two-way mixed ANOVA (nitrification for three time 
points, in brackets). 

doi:1 0.1 371/journal.pone.0086388.g004 
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Figure 5. Soil microbial community size and ammonia-oxidiz- 
ing community structure. Bars represents means ±1 SE (n = 4). 
Open bars, control treatment (NPK); grey bars, biochar treatment 
(BC3N). DNA, soil DNA content. AOA, archaeal amoA copy numbers; 
AOB, bacterial amoA copy numbers. AOA and AOB copy numbers are 
presented on soil dry mass and DNA basis. P values from two-way 
ANOVA (see Table 2). 
doi:10.1371/journal.pone.0086388.g005 

soil bacteria and archaea [72]. Soil prokaryotes typically have 
mean cell sizes of 0.6 to 0.8 |J,m in diameter [73-74] with a length: 
diameter ratio of 1.5:1 [75] and were shown to inhabit pores 
greater than 0.6 |J,m [76]. Soil fungi have slightly larger mean 
diameters, of between 1.2 to 4 |J,m [74,77] but mycorrhizal hyphae 
vary enormously in diameter ranging between 2 to 20 |J,m [78]. As 
a consequence of the high proportion of colonizable pores, 
addition of Eucalyptus charcoal to a coarse textured soil increased 
habitable pore space surface area by 1 1 % at a biochar rate of 25 t 
ha~' [70] and habitable pore volume by 6.5% at 50 t ha~' [79]. 
Indeed we found an increase of pore volume of 8% at biochar 
rates of 72 t ha~ . More importandy it has been shown that 
interior surfaces and pores are rapidly colonized, microbiota 
showing no drflFerences in densities between outer and internal 
biochar surfaces 3 years after application [79]. It has been 
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Figure 6. Relationship of archaeal (AOA) and bacterial amoA abundance (AOB) with gross soil nitrification rates. The relationships 
between archaeal amoA copy numbers (AOA) and bacterial amoA copy numbers (AOB) with gross nitrification rates were marginally significant for 
AOA (p = 0.051 ) and AOB (p = 0.072). Symbols represent the four treatments measured in September 201 2. NPK and BC3N were also measured in June 
2012. Values are means ±1 SE (n = 4). 
doi:10.1371/journal.pone.0086388.g006 



suggested that soil microbes are protected from being grazed upon 
their natural predators by colonizing pore habitats in biochar that 
are inaccessible for soil microfauna (Fig. 2). By this mechanism 
biochar pores can act as a safe microhabitat for various soil 
microbes through the exclusion of larger predators [80-81] such as 
nematodes, mites and coUembolans [73] but also microfaunal 
protozoan grazers and predators, i.e. flagellates, amoebae and 
cUiates [82-83]. Our scanning electron microscopy data suggest 
that there was an abundance of colonizable pores, particularly in 
the 0.3 to 40 ^,m range with maxima around 0.7, 3 and 15 |J.m. 
The abundant sack like structures of parenchyma cells and the 
long cave like structures of xylem vessels (Figure SI) could 
therefore offer physical protection for soil prokaryotes and fungi 
from grazers such as nematodes and amoeba, which are just 
physically too big to enter the cave or pore. This in turn could 
reduce microfaunal grazing, which leads to a decrease in 
mineralization of organic N via the soil "microbial loop" [84]. 
This scenario provides a possible explanation for the observed 
increases in nitrification rate, with bacteria, fungi, AOA and AOB 
taking refuge in the biochar. In the biochar pores AOA and AOB 
have ample access to carbon dioxide and available ammonia, 
given the high pH of this soU; this allows them to proliferate, but 
protected from grazing, thus not contributing significandy to N 
mineralization. Interestingly this scenario would suggest possible 
top down control of nitrification, with predators having been 
shown to consume up to 60% of bacterial production in some soils 
[85]. In addition to the provision of protected habitat provided by 
biochar co-location of nutrients and water on biochar surfaces and 
in pores may increase substrate availability, and thus may further 
stimulate microbial colonization [72] and nitrification. 

Conclusions and Outlook 

The main findings of this field study are that application of 
wood biochar led to marked decreases in soil organic N 
transformation rates while soil nitrification was significantly 
stimulated in this type of agricultural soil. This field study 
therefore indicates a de-coupling of the intrinsic organic and 
inorganic soil N cycles in the short- to medium-term (four to 18 
months) as a consequence of biochar addition. Assuming these 
effects persist, the deceleration of soil organic N cycling ultimately 



leads to the dichotomy of two intertwined agronomic processes; 
the build-up of soil organic N and the slowing-down of inorganic 
N release from soil organic N essential for crop growth. Addition 
of fertilizer-N in combination with biochar effectively decouples 
these processes, with plants and microbes utilizing fertilizer-N for 
growth, in turn fueUing the belowground build-up of soil organic 
N. Since soil organic N and particularly refractory soil organic 
matter are peptide dominated materials [86] , the accumulation of 
soil organic N is expected to enhance soil carbon sequestration due 
to increases in non-biochar derived soil organic matter. However, 
it remains an open question whether adsorptive capacity of 
biochar will saturate and pore-filling of biochar slows or curtails 
this process, though recent investigations with aged and fresh 
biochar suggested that the effects of biochar on organic compound 
sorption may be long lasting [87]. 

Stimulation of soil nitrification can have positive effects on plant 
N availability and crop growth but concurrently increases the risk 
of soil N losses through nitrate leaching and gaseous losses by 
denitrification and nitrifier-denitrification. The decrease in soil 
nitrate concentration and increase in microbial nitrate utilization 
however point to a relatively conservative N cycling in biochar- 
treated arable soils. It will therefore be of great importance to 
study the responses of soil organic N, nitrification, soil N losses and 
plant growth over longer time periods and in various agricultural 
settings, to better understand the long-term effects of biochar 
amendments on soil N cycling. 

Given the central role of soil organic N cycling and nitrification 
in controlling global terrestrial N cycles, the implications of this 
work are of major importance. If it is more broadly proven that 
addition of biochar combined with fertilizer-N leads to a greater 
than sum increase in non-biochar derived soil organic matter, as a 
consequence of the mechanisms described here, it may prove to be 
a fundamental factor in future soil management: as buUding soil 
organic matter through carbon farming is the key step required to 
kick start the next black/ green revolution. 

Materials and Methods 

Description of Field Trial 

The field experiment was established in Traismauer, Lower 
Austria, Austria (48° 19' 886" N, 15° 44' 308" E; 268 m a. s. 1) in 
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March 20 1 1 by the Austrian Institute of Teclmology (AIT) TuUn, 
Austria. Tlie temperate climate shows strong cominental influence 
(Pannonian climate), with mean annual precipitation of 550 mm 
and mean annual temperature being 10°C. The arable soil was a 
deep sandy to loamy silt (18.3% sand, 57.2% silt and 24.5% clay), 
classified as a calcareous Chernozem on loess. Soil pH of non- 
biochar treated soil was 7.5 in 10 mM CaCl2, the carbonate 
content was 15.8% and humus content 2.4% (Table 1). Cation 
exchange capacity of the soil was 22.5 cmol kg dominated by 
Ca^"^ (-90%) and Mg^"^ (~7%), with a small contribution by K"^ 
(2%). The site was used for agriculture for >100 years, with spring 
barley grown in 2011 and sunflower in 2012. Hardwood-derived 
biochar (80% beech, rest derived from other hardwoods) was 
produced by pyrolysis at 500°C for two hours and was purchased 
from S.C. Romchar S.R.L. (Harghita, Romania). This biochar 
had a pH of 8.6 in water and 7.5 in 10 mM CaClj and comprised 
the following elemental concentrations: 80.3% C, 9.9% O, 0.4% 
N, and 1.6%H contents (Table 4). Volatiles made up 11.1% and 
ashes 15.2% of biochar dry mass (both determined at 550°C, 
Table 4). Biochar was applied once to the field at a dry mass rate 
of 0, 24 and 72 t ha ' in March 201 1 and incorporated to a depth 
of 1 5 cm. NPK fertilizer was applied annually in spring providing 
an annual N input of 120 (2011) and 75 kg ha"' (2012), with 
corresponding annual rates of P input of 26 (2011) and 31 kg 
ha"' yr"' (2012) and K input of 50 (201 1) and 100 kg ha"' yr"' 
(2012). The four treatments were: control i.e. NPK (no biochar 
application, NPK added), BCIN (24 t biochar ha"', NPK added), 
BC3N (72 t biochar ha"', NPK added) and BC3 (72 t biochar 
ha ', no NPK added). The four treatments were arranged 
according to a Latin square and replicated four times in circular 
plots of 6.5 m diameter. 

The field is privately owned and not protected in any way. The 
local winegrower association Traisen valley (Weinbauverband 



Table 4. Element contents of beech wood biochar pyrolyzed 
for two hours at 500 "C. 




Element 


Content (% d.w.) 


C 


80.30 


0 


9.88 


Ca 


3.75 


H 


1.60 


Si 1.58 


Fe 


0.79 


K 


0.73 


N 


0.40 


Al 


0.243 


Mg 


0.198 


Mn 


0.091 


P 


0.076 


Si 0.053 


CI 


0.010 


Zn 


0.009 


Cr 


0.005 


Co 


0.002 


Ashes (550 C) 


15.2 


Volatiles (550=C) 


11.1 



Traisental) gave us the written permission to conduct the field 
experiment at the premises owned by Rudolf Hofmann who also 
agreed in his role as landowner to set up the experimental plots on 
his land. No governmental permits were required for the described 
field study, and the study did not involve endangered or protected 
species. 

Soil sampling was performed several times in 2011 and 2012 
(see Table S2 for sampling dates and for the treatments sampled 
and the pools and processes measured at each date). Soil samples 
were taken from the topsoil layer (0-10 cm) using a soil corer with 
8 cm diameter. From each treatment plot three soil samples were 
collected randomly and pooled to form one composite sample. 
The fresh samples were then weighed for bulk density calculation, 
sieved to pass 2 mm and stored at 4°C for a maximum of 2-3 days 
prior to analyses. 

Soil Physicochennical Analyses 

General soil properties. Soil water content (SWC, in g 
H2O g ' dry soil) was determined gravimetrically after drying 
fresh soil at 105°C for 24 hours in a forced-air drying oven. Bulk 
density (BD, in g cm"'^) was calculated from the fresh weight 
divided by soil core volume times dry weight: fresh weight ratio. 
The percentage of water-filled pore space (%WFPS) was 
calculated as follows: 



VoWFPS-- 



SWCxBDx 100 



BP 
2.65 



doi:l 0.1 371 /journal.pone.0086388.t004 



where 2.65 g cm"' represents the mean density of rock. 

Soil total C and N contents were measured on oven-dried 
(105°C) and homogenized soil samples. Samples were ground to a 
fine powder with a ball mill (MM2000, Retsch, Germany) and 
aliquots of 10 mg weighed into tin capsules. Total soil C and N 
contents were determined by an elemental analyzer (EA 1110; CE 
Instruments, Italy). For quantification of the soil organic carbon 
content samples were pretreated with excess acid (2 M HCl) for 
60 min and then dried at 105°C to remove carbonate from the 
calcareous soU. Other basic soil properties, such as soil pH, soil 
organic matter content, carbonate content, soil texture, total P 
content, cation-exchange capacity, base saturation, and exchange- 
able cations (see Table 1) were only measured for control (NPK) 
and biochar-treated soils (BC3N) at the Institute of Sustainable 
Plant Production, Federal Agency for Food Security (AGES) in 
Vienna, Austria, according to standard ONORM procedures 
(www.austrian-standards.at/ en). 

Biochar characterization. The biochar (beech wood char- 
coal) was characterized for element, volatiles and ash contents by 
the German Biomass Research Center (DBFZ, Leipzig, Germany) 
and results are presented in Table 4. Biochar pH was measured in 
an aqueous slurry (5% (w/v)) using a pH electrode. For 
measurements of the porosity (pores >300 nm diameter) and the 
pore size distribution of the biochar, charcoal samples were 
mounted on 0.5" aluminum specimen stubs equipped with SEM- 
Carbon foils (Agar scientific, UK). Since the specimens were 
conductive no additional coating was performed after air drying at 
40°C for 24 hours. Micrographs were taken on a Philips XL 20 
scanning electron microscope (SEM; Philips, Eindhoven, The 
Netherlands) at 350 x and 2000 x magnification with an acceler- 
ation voltage of 15 kV. For pore size and pore coverage 
measurements the software Scandium Version 5.0 (Soft Imaging 
Systems GmbH, Olympus, Vienna, Austria) and the image 
analyzing software GSA Image Analyzer (GSA, Gesellschaft fiir 
Softwareentwicklung und Analytik mbH, Germany) were used. To 
estimate pore size distribution in total, 77 15 pore measurements in 
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6 longitudinal and in 7 cross-sectional images covering 0.35 mm 
and in 25 longitudinal images covering 0.0112 mm^ were 
performed. Pores down to a diameter of 300 nm could be safely 

measured. Pore size distributions were converted to counts mm~^ 
or circular area in |Im^ mm~^ and summed up for each 1 |J.m size 
class. 

Pool sizes of extractable organic C, inorganic N and 
organic N forms. For measurements of extractable C and N 
pools fresh soils (aliquots of 4 g) were extracted for 30 min with 
30 ml 0.5 M K2SO4, subsequently filtered through ash-free 
cellulose filters and extracts kept frozen at — 20°C for later 
analyses. Total dissolved N and dissolved organic carbon 
(measured as non-purgeable organic carbon) were measured by 
high temperature catalytic oxidation on a Shimadzu TOC-VCPH 
with TNM-1 unit and ASI-V autosampler (Shimadzu Austria, 
Korneuburg, Austria). Ammonium and nitrate were determined 
with colorimetric methods, NH.^^ with a modified indophenol 
reaction and NOs" with the VCl^/Griess assay, as described in 
detail by Hood et al. [88]. Amino acid concentrations were 
quantified by a modified fluorometric OPAME procedure based 
on Jones et al. [89], optimized for free amino acid measurement in 
protein hydrolysates. The modified OPAME reagent was made up 
of 10 mg o-phthaldialdehydc dissolved in 1 ml of HPLC grade 
methanol to which 20 |J,1 of 3-mercaptopropionic acid was added. 
This reagent was mixed with 40 ml potassium borate bulfer 
(0.2 M; pH 9.5) and left to stand overnight. Sample afiquots 
(50 |J.l) were mixed with 200 |a,l reagent in the microtiter plate and 
measured after 10 min at an excitation wavelength of 340 nm and 
emission wavelength of 450 nm with a fluorescence microplate 
reader (Tecan Infinite M200, Tecan Austria GmbH, Grodig, 
Austria). As NH4"^ in the soil extracts also gives a fluorescence 
signal, NH4"'" concentration had to be determined [88] and its 
fluorescence contribution was determined via a NH4''' concentra- 
tion series and subtracted from the total fluorescence of the 
samples. Total free amino acid-N concentration was then 
calculated using glycine as a standard and assuming that each 
mol of amino acid contains 1.37 mol N. Proteins were quantified 
in K2SO4 extracts by acid hydrolysis and subsequent amino acid 
analysis (modification of [89-90]; see amino acid determination 
above). Protein hydrolysis to free amino acids was performed in 
6 M HCl and heating was performed for 15 hours at 100°C in a 
drying oven with Bovine serum albumin (BSA containing 15.8% 
N) as internal standard. Samples were then dried in an N2 stream 
and redissolved in deionized water. Determination of the NH4"'" 
concentration to correct the amino acid fluorescence was done as 
in the case of amino acid determination, with the exception that 
the sodium nitroprusside solution (the color reagent) contained 
1.5 M NaOH instead of the original 0.3 M NaOH to neutralize 
residual acid in the hydrolysates. 

Isotope Pool Dilution Assays 

Gross rates of N transformation processes were determined 
using the '^N isotope pool dilution (IPD) technique [31,91]. This 
method entails the addition of an isotopic tracer to a target pool 
and measuring the isotopic composition and concentration over 
time which allows estimating gross rates of influx to (production) 
and efflux from (consumption, immobilization) the target pool. In 
order to avoid the stimulation of consumption rates due to 
excessive addition of ' 'N labeled substrate to the respective target 
pool, each pool size was determined prior to labeling the soil 
samples. Thereafter to each target pool 10% of the pool size was 
added in the form of '"'N tracer (98 at"/) '''N). Before starting the 
IPD assays, all soils were adjusted to 40% water holding capacity 
and then the '^N tracer solutions were added to bring the soils to 



70% water holding capacity. After addition of the tracer to the 
soils the vials were shaken vigorously to homogenously distribute 
the label. The incubations were performed at room temperature 
(~20°C) and terminated by adding 0.5 M KgSO.j and subsequent 
extraction for 30 min on a rotary shaker followed by filtration 
through ash-free filter papers into new vials. Extracts were 
immediately frozen and kept for further analyses. An overview 
of the soil preparation pipeline for these isotope pool dilution 
assays is presented in Figure S2. 

Gross N mineralization and NH4^ immobilization. Each 
2 g soil sub sample was weighed into disposable HDPE vials 
(20 ml scintillation vials) and labeled witii '^NH4C1 (98 at% 
purchased from Isotec-Sigma Aldrich, Vienna, Austria). The 
incubations were terminated after 4, 24 and 48 hours by extraction 
with 15 ml 0.5 M K2SO4. Prior to isotopic analyses ammonium 
concentrations in the labeled extracts were measured photomet- 
rically. Ammonium was isolated from the extracts by microdiffu- 
sion [92]. In brief, this method implies the conversion of 
ammonium into ammonia gas which is then trapped on an 
acidified cellulose filter disc enclosed in a sc-mi-permeable Teflon 
membrane (acid trap), where the ammonia gas dissolves and 
redissociates to ammonium. The interconversion of ammonium 
and ammonia is a consequence of the pH-dependent equilibrium 
between these two N forms. To each 10 ml of soil extract one acid 
trap and 100 mg MgO were added to raise the pH to >9.5, the 
vials were closed immediately and left to stand at room 
temperature for three days. In cases where ammonium concen- 
trations exceeded the optimal range of N (10-25 (Xg NH4"*"-N 
10 ml~') for elemental analyzer-isotope ratio mass spectrometry 
(EA-IRMS) measurement, extracts were diluted with 0.5 M 
K2SO4 prior to microdiflfusion. After microdiflFusion acid traps 
were transferred to 2 ml HDPE screw-cap vials and dried in an 
evacuated desiccator over concentrated sulphuric acid for two 
days. The cellulose discs were retrieved and transferred into tin 
capsules and subjected to EA-IRMS. Conversely extracts with 
ammonium-N concentrations less than 2 |Ig NH4"''-N 10 ml~' 
were prepared to be measured with purge-and-trap isotope ratio 
mass spectrometry (PT-IRMS) [93]. For this purpose acid traps 
were made from glass fiber filters. Acid traps to be measured with 
PT-IRMS were opened and the filter discs transferred to 1.5 ml 
HDPE snap-cap vials containing 1 ml deionized water and shaken 
on a Vortex to redissolve the trapped ammonium. The alkaline 
persulfate oxidation with subsequent VClj/azide reactions were 
then appKed to convert NH4"^ via NO3 to N2O for analysis by 
PT-IRMS (see below) [93]. 

Gross nitrification and NO3 immobilization. Soil label- 
ing was done with K'^NOg (98 at% purchased from Isotec- 
Sigma Aldrich) and carried out as a mentioned above for gross N 
mineralization. As in the case of gross N mineralization a 
reduction-microdifFusion method was used to isolate nitrate-N 
after conversion by Devarda's alloy to ammonium-N from the 
extract into acid traps [92,94] and nitrate concentrations were 
determined before starting the microdiflfusions. In a first step 
ammonium was eliminated from the extracts by adding one acid 
trap and MgO (100 mg) to 10 ml sample, then the vials were 
closed and microdiffusion proceeded for three days at room 
temperature. After retrieving this acid trap, a new one was added 
together with 50 mg Devarda's alloy (45% aluminum, 50% copper 
and 5% zinc), a reducing catalyst converting nitrate to ammoni- 
um, to the soil extracts. After five days of microdiffusion at room 
temperature to collect the ammonia produced from nitrate the 
second set of acid traps was treated as mentioned above for gross 
N mineralization and isotope ratios measured via EA-IRMS (see 
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below) to estimate gross nitrification and gross nitrate immobili- 
zation. 

The organic N-IPD assays described below are novel and were 
applied here for the first time to measure gross protein and free 
amino acid fluxes in soils (but see Wanek et al. [32] for the first 
amino acid IPD by gas chromatography-mass spectrometry). 

Gross amino acid production and immobilization. The 
free amino acid pool was labeled by addition of a '^N-labeled algal 
amino acid mixture (U-'''N-98%, Cambridge Isotope Laboratories 
Europe, Radeberg, Germany). For this, each 4 g fresh soil were 
weighed into 50 ml HDPE centrifuge tubes and the '''N tracer was 
added. Samples were extracted for half an hour with 30 ml 0.5 M 
K2SO4 after 5 and 30 minutes incubation time. An ultrafiltration 
method was used to separate low-molecular weight organic N 
from high molecular weight organic N. Amicon Ultra-2 pre- 
launch centrifugal filter devices (Millipore Corporation) with a 
nominal molecular weight limit (NMWL) of 3 kDa were taken to 
collect the low molecular weight organic N fraction (mainly amino 
acids and oligopeptides) together with the inorganic N forms 
(NO3 and NH4^) in the permeate and the high molecular weight 
organic N forms (mainly proteins) in the retentate. Ultrafiltration 
was performed on a Beckman centrifuge (Beckman J2-21, 
Beckmann Coulter, Krefeld, Germany). From the labeled amino 
acid extract aliquots of 1.7 ml were centrifuged at 4000^ for one 
hour at room temperature with a swinging bucket rotor (Beckman 
JS-13.1). In order to get rid of inorganic N from low molecular 
weight organic N, each 1 ml of permeate was filled into a 2 ml 
cr}'()-vial and 20 mg Devarda's and 20 mg MgO were added 
together with an acid trap. The closed cryo-vials were kept at 37°C 
for 3 days; thereafter the acid traps were discarded. Via alkaline 
persulfate oxidation and VCls/azide reaction the low molecular 
weight organic N was then converted to N03~ and finally to N2O 
for analysis by PT-IRMS (see below). Moreover free amino acid 
concentrations were determined in each extract using the modified 
OPAME procedure as detailed above, to correct the isotope 
enrichment of free amino acids for the contributions of other low 
molecular weight organic N forms using a simple two-source 
isotope pool mixing model [95]. Free amino acid concentrations 
and the corrected isotope enrichments were used for later pool 
dilution calculations. 

Gross protein production and consumption. Soil labeling 
was performed as in the case of the amino acid-IPD assays but 
with the '^N tracer being a purified algal crude protein (98 
at%'^N purchased from Isotec-Sigma Aldrich) and with incuba- 
tion times of 5 min, 30 min, 4 and 24 hours. After ultrafiltration 
(similar to amino acid-IPD but centrifugation parameters altered) 
of each 1.7 ml of the labeled extract at 7500 g for one hour at 
room temperature with a fixed angle rotor (Beckman JA-20), the 
permeate was put aside. The high molecular weight organic N 
containing retentate was washed twice with each 1 ml 50 mM 
K2SO4 via 30 min centrifugation at room temperature and 
7500 g, to remove inorganic N and low molecular weight organic 
N. Thereafter the retentate was recovered by 2 min centrifugation 
at 1000 g at room temperature and diluted to 0.5 ml with 50 mM 
K2SO4. The alkaline persulfate oxidation and subsequent VCI3/ 
azide reactions were then used to convert the high molecular 
weight organic N via NO3 to N2O for analysis by PT-IRMS (see 
below). 

Alkaline persulfate oxidation (APO) of N„rg (NH4*) to 

NO3 . The APO reaction was used to convert ammonium and 
organic N forms into nitrate [93]. At elevated temperature 

persulfate decomposes forming persulfate radicals that efficientiy 
oxidize organic and inorganic compounds, in the case of organic N 
and ammonium producing nitrate. The persulfate reagent was 



prepared by dissolving 44 g Na2S208, 16.8 g NaOH and 30 g 
H3BO3 in deionized water with an end volume of 1 L [96]. For 
digestion aliquots of 750 |Xl of sample were pipetted into 1.5 ml 
HPLC glass vials and each 750 |J,1 of persulfate solution were 
added. Vials were closed immediately (to avoid losses of 
ammonium due to the alkaline pEI) and autoclaved for one hour 
at 120°C. Digested samples were stored at 4°C until further 
analyses. 

VCls/azide conversion of NOs^ to N2O. The isotope ratios 
of nitrate after persulfate digestion were measured by conversion 
of nitrate to N2O by a coupled VCI3/ azide reaction in which in a 
first step, VCI3 is used to reduce nitrate to nitrite, which in a 
second step is further reduced to N2O by sodium azide [93]. For 
this, aliquots of 1 ml sample were pipetted into 12-ml screw cap 
exetainers with butyl rubber septum (IVA-VC309; IVA Analyse- 
technik, Merbusch, Germany) which were closed tightly and 
purged with hehum for 10 min to get rid of air-N20 in the 
headspace. Thereafter 160 |J,1 sodium azide and finally 1 ml VCI3 
were added using gas-tight syringes by injection through the septa. 
Exetainers were put on a rotary shaker and kept at 37°C for 18 
hours. The reaction was stopped by neutralizing and deactivating 
the sodium azide through injecting each 160 (il 6 M NaOH into 
the exetainers and N2O was analyzed by PT-IRMS. 

PT-IRMS. The exetainer vials were then loaded into a 96-slot 
autosampler with a double-hole needle (GC-PAL, CTC Analytics, 
Zwingen, Switzerland) connected via a Gasbench II lu-adspace 
analyzer (Thermo Fisher, Bremen, Germany) and a PreCon 
cryotrapping device (Thermo Fisher) to a gas isotope ratio mass 
spectrometer (IRMS, Delta V Advantage, Thermo Fisher) [93]. 
Through the double needle of the autosampler an auxiliary helium 
flow purged the contents of the exetainer vials through a water and 
CO2 scrubber into a sequence of two cryo-traps held at liquid N2 
temperature, and the N2O was then released and residual CO2 
and air (O2, N2) separated by gas chromatography. The purified 
sample N2O was then transferred into the mass spectrometer via 
an open split with repeated pulses of standard N2O gas for internal 
calibration (Air Liquide, Vienna, Austria). For calibration of 
at% N respective standards varying in 'N: N ratios were run 
throughout the sample preparation methods. For further details on 
PT-IRMS analyses see Lachouani et al. [93]. 

EA-IRMS. Acid traps and homogenized dry soils were 
transferred into tin capsules and directiy analyzed for N content 
and at%'^N by EA-IRMS. The system consisted of an elemental 
analyzer (FA 1110, CE Instruments, Milan, Italy) connected via a 
ConFlo III interface (Thermo Fisher) to the IRMS (Delta'"^"'^, 
Thermo Fisher). External calibration was performed using 
laboratory standards for isotope composition and N concentration 
[93]. 

Isotope pool dilution calculations. Gross production (GP) 
and gross consumption (GC) rates, i.e. influx into and efflux from 
the target pools that was labeled with '^N in the course of IPD 
experiments was then calculated based on the isotope pool dilution 
equations developed by others [97]: 
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h-h 



In 



1 + - 



APE,^ 



\ 



In 



where C,i, C,2, APE,i and APEg are the concentrations of N (|a,g 
N g ' soil dry weight) and atom percent excess of ' 'N (APE, in %) 
in the respective target N pool at ti and tj (days), i.e. the times of 
IPD sample termination. APE was calculated as the difference 
between the at%'^N of the given N pool at time t and the 
background '''N abundance (at%'''N) measured in separate non- 
labeled samples. Gross rates are therefore expressed in |lg N g ' 
soil dry weight d~'. 

Extraction of nucleic acids and quantitative PGR 
analyses. Nucleic acids were extracted from soil samples 
(0.5 g fresh weight) using a modification of the method described 
by Griffiths et al. 2000 [98] . Samples were loaded in 2 mL Lysing 
Matrix E tubes (MP Biomedicals, lUkirch, France) and cells were 
disrupted by bead-beating with 0.5 mL of an SDS-based lysis 
buffer (1% SDS, 0.7 M NaCl, 0.1 M Na2(S03), 0.05 M EDTA 
pH8 and 0.1 M Tris/HCl pH 7.5) and 0.5 mL of a phenol:- 
chloroform:isoamyl-alcohol (25:24:1) mixture (pH 7.5). Precipita- 
tion of nucleic acids was performed for 1 h at room temperature 
with 20 |J,g glycogen (MBI Fermentas, Vienna, Austria), 0.5 mL 
isopropanol and 40 |i,L 5 M NaCl. Precipitated nucleic acids were 
washed twice in 500 |t,L ice-cold 70% ethanol, re-suspended in 
50 |XL DEPC-treated water and purified by filtration through 
PVPP mini-columns centrifuged at lOOOx^ for 2 min at 10°C. 
Purity and concentration of purified nucleic acids was checked 
spectrophotometrically using a NanoDrop spectrophotometer 
(Thermo scientific) and their integrity was determined by agarose 
gel electrophoresis. A SybrGreen I quantitative PGR approach 
was used to quantify thaumarchaeal and bacterial amoA gene 
abundance. Thaumarchaeal amoA genes were amplified with PGR 
primers Arch-amoA-104F/Arch-amoA-616R [99] and bacterial 
amoA genes with primers amoA-lF/amoA-2R [100]. Duplicate 
dilution series (10 -10 copies) of PGR products obtained with 
M13F and M13R vector primers from Topo® TA cloning® 
vectors (Invitrogen) containing i;ither a 632 bp fragment of 
Mitrososphaera viennensis EN76 [101] amoA gene, obtained with 
PGR primers Arch-amoA-7F/Arch-amoA-638R [99], or a 491 bp 
fragment of Nitrosospira multiformis ATGG 25196 amoAl gene 
(Nmul_A0799/Nmul_A2325), obtained with primers amoA-lF/ 
amoA-2R, were used as standards for the thaumarchaeal and 
bacterial amoA assays, respectively. Thaumarchaeal amoA standards 
were amplified with an efficiency of 93% and / value of 0.997 and 
bacterial amoA standards with an efficiency of 93% and / value of 
0.991. All reactions (20 |iL) contained 10 nL of 2x Quantifast™ 
SYBR® Green PGR Master Mix (Qiagen, Vienna, Austria), 
0.2 mg.mL"' BSA, 1 (Archaea) or 0.8 (Bacteria) |lM of each 
primer and for the analysis of soil samples 10 ng of purified nucleic 
acid extracts. Amplification conditions were 95 °G for 15 min, 
followed by 40 cycles of 95°C for 15 s, 60°G for 30 s (combined 
annealing and cxtc-nsion step) and plate read at 78°G for 
thaumarchaeal amoA genes and 35 cycles of 95°G for 15 s, 62°G 
for 30 s and plate read at 82°G for bacterial amoA genes, before a 
final extension step at 60°G for 10 min. Reactions were carried out 
in a Mastercycler® RealPlex^ thermocycler (Eppendorf Austria 
GmbH, Vienna, Austria) and the specificity of the amplification 



was assessed by melting curve analysis and agarose gel electro- 
phoresis. 

Statistical analyses. Mixed between-within subjects ANO- 
VAs (split-plot ANOVA) followed by Tukey HSD post hoc tests 
were performed in SPSS version 16.0 for Windows (SPSS Inc., 
USA) to test for the effects of treatment and sampling time on 
biological and physicochemical soil parameters. If normality 
(Shapiro-Wilk test) and variance homogeneity (Levene's test) were 
violated, data were log, srjuare root or reciprocal transformed to 
meet these assumptions. Mauchly's test was used to test the 
assumption of sphericity and if significant, the degrees of freedom 
of the F-distribution were corrected using the Greenhouse-Geisser 
procedure as the estimated epsUon was always less than 0.75. 
Biochar effects on pools and processes that were only measured 
once were tested by oneway ANOVA. Regressions and correla- 
tions were calculated with Statgraphics Genturion XVI (StatPoint 
Technologies Inc., USA). Simple least squares linear regressions 
were calculated to estimate the amount dependency of biochar 
effects, with added biochar amount as independent variable. 
Gorrelations between microbial community structure parameters, 
or with nitrification rates were performed based on the Pearson 
product-moment correlation procedure. 

Supporting Information 

Figure SI Structure of beech wood biochar shown by 
SEM analysis. Gross section across early to late wood transect in 
beech wood charcoal (top left), longitudinal section through beech 
charcoal showing open xylem vessels and xylem parenchyma cells 
(top right), close-up of xylem pits (bottom left) and close-up of a 
plasmodesmata channel (bottom right). 
(DOGX) 

Figure S2 Scheme of the workflow for sofl sample 
preparations to measure organic and inorganic N 
transformation rates using ^^N isotope pool dflution 

approaches. 

(DOGX) 

Table SI Results of ANOVA of soil physicochemical 
factors and gross nitrification rates for the factors 
treatment and sampling date for all four treatments. 

(DOGX) 

Table S2 Sofl sampling dates, treatments sampled and 
measurements done on these samples. 

(DOGX) 
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